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DISTRIBUTION  LIST 


ABSTRACT 


ANTI  is  a  computer  program  which  scales  existing  neutron  transport  data 
given  at  sea  level  to  determine  the  energy  spectrum  of  neutrons  transported 
fTom  a  given  source  to  a  receptor.  The  energy  spectrum  can  be  determined 
over  wide  ranges  of  altitude  and  distances  of  separation  between  source 
and  receptor#  The  basic  transport  data  was  obtained  from  Moments  Methods 
runs  by  Kaman  Sciences  Corporation  and  Monte  Carlo  results  reported  by 
Braddoc'k,  Dunn  and  McDonald. 
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INTRODUCTION 


The  computer  program  ANTI  was  written  at  the  request,  of  the  Nuclear 
Effects  Branch1  ,  Engineering  Sciences  Laboratory,  (ESL),  Feltman  Research 
Laboratory.  The  model  was  formulated  by  Dr.  F.  J.  Jankowski  who  was  a 
consultant  to  ESL. 

The  code  was  developed  to  provide  an  inexpensive  but  reliable  means 
of  obtaining  the  flux  distribution  caused  by  a  neutron  source  in  the  atmos¬ 
phere.  ANTI  results  were  compared  against  the  original  source  of  the 
transport  data  (ref  1),  as  well  as  against  various  other  air  transport 
studies.  These  checks  gave  quite  acceptable  agreement, 

ANTI  computes  the  flux  and  energy  spectrum  of  neutrons  transported 
from  a  given  source  to  a  receptor.  The  difficult  problem  of  neutron  trans¬ 
port  is  handled  by  scaling  existing  transport  data  given  at  sea  level  (see 
ref  1).  The  transport  table  provides  the  code  with  previously  calculated 
flux  data  out  to  a  distance  of  6500  feet  at  sea  level  (approx.  17  mean  free 
paths  for  a  lii  Mev  neutron).  The  flux  data  is  tabulated  for  33  discrete 
ranges  from  0  to  6500  feet.  Interpolation  techniques  make  it  possible  to 
find  the  flux  at  any  given  range  less  than  6500  feet.  For  longer  distances 
an  extrapolation  procedure  is  employed. 

The  program  can  list  in  a  convenient  format  the  complete  set  of 
transport  data  built  into  the  program.  The  transport  table  may  be  expanded 
or  replaced  to  include  additional  data  should  it  become  available. 

ANTI  is  written  in  FORTRAN  IV  for  the  IBM  360/65  and  the  CDC  6500 
computers. 

Nuclear  Effects  Br.  is  now  the  Threat  and  Susceptibility  Br*,  Concepts 
and  Effectiveness  Division,  Nuclear  Development  and  Engineering  Directorate 
(NDED), 
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MATHEMATICAL  MODEL 


The  computer  program  ANTI  is  used  to  calculate  the  flux  and  energy- 
spectrum  of  neutrons  which  a  target  (receptor)  receives  from  a  given 
nuclear  detonation  (source).  The  geometry  is  illustrated  in  Figure  1, 


FIGURE  1  (Geometry  of  Problem) 

where 

hfc  is  the  altitude  of  the  source, 

hr  is  the  altitude  of  the  receptor. 

D  is  the  distance  between  the  source  and  receptor. 

<P  is  the  angle  between  the  vertical  axis  and  the  vector  from  the 
source  to  the  receptor. 

The  source  spectrum  is  specified  for  a  sat  of  discrete  energy  levels 
{SL(1),  i  =  l,n}.  The  number  of  neutrons  {Ng(i)}  is  specified  at  each 
energy  level  or  within  each  interval  of  the  given  spectrum. 

As  the  neutrons  travel  from  the  source  to  the  receptor,  they  inter¬ 
act  with  the  atmosphere;  this  interaction  is  a  function  of  the  distance 
traveled  and  the  density  of  air  along  the  trajectory  of  the  neutrons.  The 
program  assumes  the  attenuation  along  the  given  path  (not  including  the  in¬ 
verse  square  attenuation)  is  the  same  as  the  attenuation  the  neutrons  encounter 
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in  traveling  a  distance  D0  at  sea  level.  From  reference  (l)  D0  is  computed  as 

D 


(?(r)dr 


(1)  "o 

where 

is  the  density  of  air  at  sea  level. 

is  the  density  of  air  at  the  distance  T  from  the  source. 
In  ANTI,  the  density  of  air  is  a  function  of  altitude  only. 
(1)  is  computed  as 


(2) 


*P(h)dh 

<hb 


using  the  transformation 

r  -  — h— 

(3)  COS0 

9 

The  density  of  air  is  built  into  the  program  for  discrete  values  of 
altitude  which  range  from  0  ft  to  200,000  ft  in  increments  of  10,000  ft 
(reference  2).  Between  two  adjacent  altitudes,  the  density  is  assumed  to 
be  of  the  form 

(h)  Pw  -&V  EXP  {-^-N  | 

where 


(°(h) 

i 


is  the  density  of  air  at  the  altitude  h. 
denotes  the  interval  in  the  table  of  air  densities  in  which  h  lies. 


k 


hi  is  the  lower  altitude  of  the  i"^  interval. 

Si  is  inverse  of  the  slope  of  log  (h)  in  the  i^h  interval. 

In  general,  the  source  and  receptor  may  be  located  in  different  intervals. 
Therefore,  it  is  necessary  to  carry  out  the  integration  in  (3)  as  the  sum  of 
integrals  over  each  interval.  That  is 

,  phi+l  ^  phj+l  nhr 

(5)  d  -  ±—  p(h)dh  ♦  >  P(h)dh  +  ^(h)dh  ! 

°  j4lh3  \  j 

where  NK  is  the  interval  in  which  the  receptor  lies  and  the  source  is  in 
the  lkh  interval. 

Within  the  jth  interval,  equation  (li)  can  b*-  integrated  analytically  as 


It  is  obvious  that  in  an  actual  case,  the  only  variable  in  equation  (6) 
which  changes  if  the  integration  is  over  a  full  interval  is  ft  Therefore, 

COS  D  ( j )  was  comnu'iv’d  for  each  interval  in  the  air  density  table  and  built 

o 

into  the  program  as  a  table.  For  a  given  engagement,  ANTI  first  determines  which 

intervals  the  source  and  receptor  are  located  within  and  evaluates  the  first 

and  last  integrals  in  (lJJ.  The  integrals  in  the  sum  are  obtained  by  multiplying 

the  tabular  values  of  D0(j)  cos^  by  l  • 

COS^b 

Transport  data  from  reference  1  and  2  are  built  into  the  program 
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in  the  form  of  tables.  The  tables  givo  the  probability  of  neutrons  (flux 
multiplied  by  iitrD^3)  .being  transported  from  the  source  to  the  receptor  as 
a  function  of  source  energy  ST(j),  receptor  energy  interval  RT(i),  and 
distance  between  the  source  and  receptor  R.  The  values  for  R  ■  10li2  m 
are  illustrated  in  Table  1.  Data  is  given  for  33  values  for  R  between  0 
ft  and  6500  ft.  The  distances  were  selected  so  that  the  spacing  was  fine 
close  to  zero  and  increased  at  the  distances  close  to  6500  ft.  The  reason 
for  this  is  that  the  Dq  will  not  necessarily  be  one  of  the  tabulated  R  values 
and  logarithmic  interpolation  is  used  to  determine  the  transport  at  Bq,  if 
Dq  is  less  than  6500  ft.  If  Dq  is  greater  than  6500  feet,  logarithmic 
extrapolation  is  used  to  determine  the  transport  table  at  DQ, 

Once  Dq  and  the  corresponding  table  of  receptor  flux  (function  of 
source  and  receptor  energy  interval)  are  calculated,  it  remains  to  determine 
the  transport  of  the  given  source.  In  general,  the  given  source  spsctrum 
is  not  the  same  as  tlie  source  spectrum  in  the  transport  table.  This  is 
overcome  by  the  following  interpolation  scheme.  Let 

(SL(i)  i  »  l,n}  be  the  energy  levels  in  the  given  source  spectrum 
(or  the  upper  bound  of  the  energy  level,  for  continuous  spectra), 

{ ST( J )  j  ■  l,m}  be  the  source  energy  levels  in  the  transport  table 
(or  the  upper  bound  of  the  energy  level  for  continuous  spectra), 

Hg(i)  be  the  number  of  neutrons  within  the  i*1*1  energy  level. 

*ST(j)  is  the  upper  bound  of  the  interval  if  continuous  spectrum  is  specified, 

•«RT(i)  is  the  upper  bound  of  the  energy  interval  in  which  the  neutron 
appears  at  the  receptor. 
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The  following  conventions  are  required  in  the  program: 

SL(l)  >  SL(2)  >  .  .  .  >  SL(n) 

ST(l)  >  ST(2)  >  .  .  .  >  ST(m) 

SL(l)  =  ST(l)  (=11*  Mev  in  the  initial  program) 

SL(2)  =  SP(2)  (=12  Mev  in  the  initial  program) 

SL(n)  >  ST(m) 

The  maximum  and  minimum  energy  bounds  that  a  neutron  may  have  are 
1U  Mev  and  .001:  Mev  respectively.  That  is: 

.001:  <  SL(i)  <  Hi  MEV 

Neutrons  with  energies  below  ,00ii  Ibv  are  ignored. 

From  the  given  source  data,  the  program  generates  a  curve  which  gives 
the  number  of  neutrons  with  energy  greater  than  a  given  energy  level  (NGS(i)} 


as  a  function  of  energy  level.  From  the  restrictions  above,  the  number 

greater  than  lh  Mev  must  be  0.  NGS(i)  is  computed  as 

i 


(7) 


N0S(i) 


I  V" 


K“1 


In  equation  (7)  i,the  index  of  NGS,  refers  to  the  lower  bound  r,f  the  1^ 


energy  interval. 

The  number  of  neutrons  greater  than  the  energy  level  in  the  trans¬ 
port  table  NGT(j)  is  obtained  by  the  following  interpolation  formula 


(8) 


NGT(j)  «  NGS(i-l)  + 
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,  I 

if 

SL(i-l)  >  ST(J)  >  SL(i)  i  =  2,3, ...,n 

j  3  ?,3,...,m 

t 

The  numb or  and  intensity  of  neutro;is  at  the  source  in  terns  of  the 
source  spectrum  in  the  transport  tables  is  calculated  as 

NgjCj)  3  NQT(J)  -  NOT(j-l)  3 -  2,3, ...,m 

(9) 

NST(l)  =  NGT(l)  ■  NGS(l)  »  Ns(l) 

Equations  (7),  (8),  and  (9)  are  usod  to  convert  the  given  source  data  to 

I 

correspond  to  the  source  spectrum  in  the  transport  table.  The  neutron 
flux  received  by  the  target  from  the  given  source  is  then  calculatod  as 

(10 )  n-1 

0T(j)  °  r-ir-  I  NST(K)0(K,j  ) 
hndM  K-l 

i  i 

0(k,j )  is  the  previously  mentioned  probability  of  a  neutron 
starting  out  in  energy  interval  K  ending  up  in  energy  interval  j  after 
traversing  an  equivalent  sea  level  distance  DQ. 

D  is  the  distance  R  between  burst  point  and  receptor, 

0^(j)  is  the  neutron  flux  at  the  receptor  with  energies  in  the 
interval.  The  energy  intervals  correspond  to  the  receptor  spectrum  in 

the  transport  table.  Equations  (7),  (8),  and  (9)  can  be  used  to  convert 

this  data  to  any  desired  spectrum  for  output  purposes. 
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SUMMARY  OF  FROGRAM  CAPABILITY 


(1 )  The  source  spectrum  may  have  up  to  25  energy  intervals  (ranging 
from  lij  .ev  to  ,00lt  Mev 

(2)  For  each  energy  interval,  the  program  will  accept  either  the 
number  of  neutrons,  the  number  of  neutrons  per  kiloton  or  the  fraction 
of  the  total  number  of  neutrons  at  the  source  having  the  given  energy, 

(3)  If  this  data  is  normalized  to  one  kiloton,  then  the  number  of 
kilotons  at  the  source  must  be  specified, 

(U)  The  altitude  of  the  source  (h. ),  the  altitude  of  the  receptor 
(l^)  the  distance  between  the  source  and  the  receptor  (D)  and  the  number 
of  kilotons  (nkt)  at  the  source  can  be  parameterized  respectively  while 
reading  the  source  spectrum  once.  The  maximum  number  of  variations  are 
20  values  of  hb,  20  values  of  hr,  30  values  of  D*  and  10  values  of  nkt, 

(5)  The  maximum  number  of  cases  handled  can  be  easily  increased 
by  increasing  the  dimension  of  the  appropriate  variable  array, 

(6)  The  transport  data  at  sea  level  is  built  into  the  program  and 
can  be  modified  in  a  straight-forward  manner  as  more  accurate  data  is 
obtained, 

(7)  An  option  exists  in  the  program  to  print  the  transport  data  at 
sea  level  which  is  built  into  the  program, 

(8)  An  option  exists  in  the  program  to  print  the  transport  data  at 
sea  level  for  the  given  source  and  receptor  conditions, 

(9)  The  given  source  spectrum  does  not  have  to  be  the  same  as  the 
source  spectrum  in  the  transport  table, 

(10)  The  spectrum,  printed  out,  does  not  have  to  be  the  same  as  t he 
given  source  spectrum  in  the  transport  table. 
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DESCRIPTION  OF  INRJT 


This  section  describes  the  input  procedure  necessary  to  use  ANTI, 

The  column  headed  "Card  Set"  refers  to  one  or  more  cards  of  the  same 
type  and  the  cards  are  arranged  in  the  order  of  card  sets.  The  column 
headed  "Symbol"  gives  the  variable  name  used  in  the  program.  The 
column  headed  "Card  Column"  gives  the  column(s)  in  which  the  numerical 
value  is  punched.  The  column  headed  "Definition"  describes  the  meaning 
of  the  symbol  in  the  program.  The  column  headed  "Units"  defines  the 
units  that  the  program  assumes  for  the  variable.  Some  symbols  refer 
to  an  array  name.  In  this  instance  a  subscript  is  printed  after  the 
symbol  and  the  "Conversion"  Column  describes  the  field  width  of  the 
card.  If  more  than  one  card  is  necessary  to  describe  the  array,  additional 
cards  are  added  using  the  same  format  as  the  first  card. 
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CARD 

CARD 

SET 

SYMBOL 

COL. 

CONVERSION 

DEFINITION 

UNITS 

1 

JTABS 

1 

11 

Control  to  print  transport 
data  JTABS=  0  do  not  print 

JTABS*  1  print 

1 

NTRT 

2 

.11 

Control  to  print  transport  data 
for  specified  equivalent  sea  level 
ranges  NTRT-  0  do  not  print 

NTRT*  1  print 

1 

NRANGE 

3 

n 

Option  to  determine  range  for  a 
specified  flux  value 

NRANGE=  0  do  not  use  option 

NRANGE*  1  use  option 

2 

PHI 

1-10 

E10.3 

Only  use  if  NRANGE=1. 

Value  of  specified  1-Mev 
equivalent  flux. 

3 

NSOR 

1-3 

1.3 

Number  of  source  altitudes  to  be 
considered  (maximum  of  20 ) 

3 

NREC 

U-6 

13 

Number  of  receptor  altitudes  to  be 
considered  per  source  altitude 
(maximum  of  20) 

3 

NDSR 

7-9 

13 

Number  of  distances  between  each 
source  and  each  receptor  (Maximum 
of  30) 

EXAMPLE:  If  NS0R=2,  : 

NREC=3 ,  NDSR=ii 

,  then  2h  separate  cases  will  be  calculated. 

3 

NSIT 

10-12 

13 

Number  of  source  energy  levels 
(maximum  of  25) 

• 

3 

NSST 

13-15 

13 

Control  to  specify  source  energy 
levels  different  than  that  of  trans¬ 
port  data 

NSST*  0  energy  levels  are  the  same 
NSST  =  as  transport  data. 

1  energy  le  \els  are  different 
than  transport  data. 
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CARD  CARD 

SET  SYMBOL  COL.  CONVERSION  DEFINITION 


3  NOPD  16-18  13 


3  NNTSD  19-21  13 


Control  to  print  energy  levels 
at  the  receptor  different  than 
transport  data 

NOPD*  0  Output  spectrum  is  the 
same  as  the  transport 
data. 

N0PD=  1  output  spectrum  is  dif¬ 
ferent  than  the  trans¬ 
port  data. 

If  source  is  normalised  to  1KT, 
NNTSD  is  the  number  of  sources 
to  consider  with  number  of  KT#.. 


UNITS 


EXAMPLE:  If  the  source  is  3KT,  then  the  energy  is  determined  by 
3x( results  for  1KT) . 


3  NPCT  22-2U  13 

h  AHS0RC(I)  1-70  7F10.3 
5  SSPCTR(I)  1-70  7F10.3 


Option  to  specify  the  fraction 
of  the  total  neutron  flux. 

NPCT*  0  source  spectrum  is  specified 
as  the  neutrons  per  energy 
level . 

NPCT=  1  source  spectrum  is  specified 
as  the  fraction  of  the  total 
neutron  flux. 

The  altitudes  of  the  source  to  be 
considered  (i-l,  NSOR)  ft. 

The  energy  levels  of  the  given 

source  in  decreasing  order  (l=l,NSIT).  Mev 


NOTE:  Omit  card  sets  6  and  7  if  N0PD=0 


NOSS 

1-3 

13 

The  number  of  energy  intervals  in 
the  output  spectrum  plus  1. 

- 

OSRID(I) 

1-70 

7F10.3 

The  energy  levels  of  the  output 
spectrum  (l*l,  NOSS) 

Mev 

13 


CARD 

SET 

SYMBOL 

CARD 

COL. 

CONVERSION 

DEFINITION 

UNIT’S 

NOTE: 

Card  sets  8  and  9  have  two 

options. 

Option 

Is  If 

NPCT=1 

8 

TNNTS 

1-10 

F10. 3 

The  total  number  of  neutrons  at 
the  source 

- 

9  PCTNT(I) 

1-70 

7F10.3 

The  fraction  of  the  total  flux 
in  the  Ith  energy  interval  at 
the  source  (l=l,NSIT-l) 

Option 

2:  If 

NPCT=0 

8  ANNTS( I) 

1-70 

7F10.3 

The  neutron  flux  at  the  source 
having  energy  in  the  Ith  interval 
(1=1,  NSIT-1 ) 

Card  set  9  is 

omitted. 

NOTE:  Card  set  10  used  only  if  NNTSD  >  0 
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SSNKT(l) 

1-70 

7F10.3 

The  number  of  kilotons  for  the 

Ith  source  ( 1=1,  NNTSD) 

Kt 

11 

AHREC(l) 

1-70 

7F10.3 

The  altitude  of  the  Ith  recep¬ 
tor  (l»l,  nrk) 

ft 

T  0 

•4.(1 

DRSR(l) 

1-70 

7F10.3 

The  Ith  distance  between  the 
source  and  receptor  ( 1=1,  NREC) 

cn 

NOTE:  Card  sets  11  and  12  are  repeated  for  each  source  altitude  to  be 
considered  (i.e.  repeat  NSOR  times), 

13  RATS  1-ii  All  ’TtATS”  is  punched  in  columns  1  through 

li  and  signifies  the  end  of  all  input 
data 


111 


Caro  field 


DESCRIPTION  OF  OUTPUT 


Basically  ANTI  computes  the  flux  and  energy  spectrum  of  neutrons 
which  a  receptor  receives  from  a  given  source.  In  addition  to  this  data, 
ANTI  also  prints  a  listing  of  all  the  input  parameters,  the  complete 
transport  data  as  an  option,  and  the  transport  table  for  the  equivalent  sea 
level  range  (D0). 

The  initial  output  is  a  complete  listing  of  all  the  input  data 
(card  images).  This  is  done  by  first  storing  an  input  card  onto  a  disk 
and  then  printing  the  card  before  a  new  card  is  read.  Then  the  data  for 
an  individual  case  are  read  from  the  disk. 

If  the  option  to  print  the  complete  transport  table  is  selected, 
this  data  will  appear  next  ( JTABS=l) .  The  transport  data  for  one  value 
of  range  is  printed  on  each  page.  The  first  line  contains  the  range  value 
in  feet.  The  next  line  contains  the  source  of  the  receptor  spectrum.  The 
entries  in  the  table  are  then  the  fraction  of  neutrons  (flux  multiplied  by 
iirrD02)  transported  from  the  source  to  the  receptor.  This  data  is  repeated 
for  each  range  value  (33  pages  of  output). 

A  second  option  (NTRT=l)  is  to  print  the  transport  table  for  the 
equivalent  sea  level  distance  (D0).  The  format  for  this  output  is  the 
same  as  above  except  J^0  is  printed  on  the  first  line  of  the  page. 

Finally  the  basic  output  is  printed.  This  consists  of  the  neutron 
spectra  for  each  set  of  source  (burst)  altitude,  receptor  altitude,  source- 
to-ieceptor  slant  range  and  yield  (KT). 
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A  row  of  surmary  information  is  listed.  This  includes  values  for  each 
of  the  parameters  above  and  the  equivalent  sea  level  range  (Dn).  Next  tables 
of  flux  (nt/cm2)  and  total  neutrons  by  energy  group  are  printed.  These 
are  supplemented  by  tables  of  direct  flux  (if  NSSf  equals  zero)  or  scattered 
flux  (if  NSST  equals  zero  and  output  energy  boundaries  correspond  to  all 
twelve  source  energy  boundaries). 

NOTE:  NET  FLUX  =  DIRECT  +  SCATTERED 

The  final  output  consists  of  summations  over  all  energy  groups  of  total 
fluence,  one-Mev  equivalent  flux,  and  RADS. 


XT? 
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APPENDIX  A 


BRIEF  DESCRIPTION  OF  THE  TRANSPORT  TABLES 

The  air  transport  tables  are  built  into  the  program  as  data  statements. 
The  data  is  normalized  to  one  neutron  per  source  group.  They  are  very 
dependent  on  the  13  flux  levels  at  the  source,  17  flux  levels  at  the  receptor 
and  33  values  of  range.  If  any  of  these  numbers  are  changed,  the  dimension 
statements,  equivalence  statements,  and  data  statements  which  describe  the 
new  air  transport  table  must  be  consistent  with  the  new  numbers. 

Each  data  statement  gives  the  transport  data  as  a  function  of  range 
for  a  given  source  energy  group  and  receptor  energy  group.  The  ranges  are 
defined  in  increasing  order  and  the  energy  intervals  are  in  decreasing  order. 
The  data  statements  are  ordered  such  that  the  receptor  energy  varies  for  each 
source  energy  level. 

These  data  statements  were  generated  by  a  separate  program  which  is 
available  to  generate  an  entirely  different  table. 
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APPENDIX  B 


I 


PROORAM  SYMBOLS  AND  UNITS 


Symbol 

Definition 

Units 

ADENS(I) 

Altitudes  at  which  density  of  air  has 
been  tabulated  in  the  program. 

(21  values  flrom  0  to  200,000  ft), 

ft 

AHREC(I) 

Altitude  of  receptor.  (Maximum  of  20 
values) 

ft 

AHS0RC(  I) 

Altitude  of  source.  (Maximum  of  20 
values) 

ft 

ANNTR 

Neutron  flux  at  receptor  for  a  given 
energy  group.  (Partial  summation) 

nt/ft* 

ANNTS(l) 

Number  of  neutrons  in  the  Ith  energy- 
group  of  source  spectrum. 

BNNTS( I) 

XU 

Number  of  neutrons  in  the  I  energy- 
group  in  terms  of  source  spectrum  of 
the  transport  table. 

CSI 

Cos  ^  (see  figure  l)  >• 

DBRS 

Distance  between  source  and  receptor 
being  evaluated 

ft,  cm 

1 

DENS(l) 

Density  of  air  at  altitude  ADENS(l) 

lb/ft* 

DIRECT(l; 

XL 

Unscattered  flux  in  I  energy  group 

nt/cm* 

DNR(l) 

Sf(h)  dh  from  h  -  ADENS(l)  to 

J  h  ■  ADENS(l+l),  where S  is  the  den¬ 
sity  of  air  (20  values) 

lb/ft a 

drsr(i) 

Distance  between  source  and  receptor 
(maximum  of  30  values) 

cm 

EQUIV(l) 

One-Msv  Equivalent  weighting  factors 
(for  each  of  16  output  energy  groups  of 
the  transport  table  spectrum) 

EQUR 

Equivalent  ERSR(l)  distance  at  sea  level 

3 

O 

a 

FLUXR( I) 

Rads  equivalent  of  flux  in  the 

Rads 

output  energy*  youp. 
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Symbol 

FLUXI(I) 

FREE 

ha(i) 

FTCMSQ 

HREC 

HSCRC 

IHRI 

msi 

NDNS 

NDSl 

NNTSD 

NOFD 

NOSS 

NR 

NRANGE 

NREC 

NRLS 

NRSI 


Definition  Units 

One-Mev  equivalent  of  flux  in  the 
i  output  energy  group 

liirD*,  where  D  is  the  distance  from  ft* 

source  to  receptor 

Slope  of  curve  of  log(f  )  vs.  altitude  1/ft 

from  h  -  ADENS(l)  to  h  ■  ADENS(l+l) 

Conversion  from  ft*  to  cm*  cm*/ft* 

Altitude  of  receptor  being  evaluated  ft 

Altitude  of  source  being  evaluated  ft 


Receptor  is  located  between  ADENS(lriRl) 
and  ADENS( IHRI+1 ) 

Source  is  located  between  ADENS(lHSl) 
and  ADENS(  IHSI+l) 

Number  of  altitudes  at  which  the  density 
of  air  is  specified.  (21  altitudes  from 
0  to  200,000  ft) 

Number  of  distances  between  source  and 
receptor  to  be  investigated. 

Option  to  specify  the  number  of  kilotons 
in  the  source. 

Option  to  specify  the  output  energy  spectrum 

Number  of  energy  intervals  in  the  output 
spectrum. 

Number  of  range  points  at  which  the  transport 
table  flux  is  mpecified.  (33  points  from 
0  to  6500  ft. 

Option  to  compare  flux  at  a  distance  to  a 
specified  one-Mav  equivalent  flux. 

Number  of  receptor  altitudes. 

Number  of  energy  levels  in  the  output  spec¬ 
trum  of  the  transport  table  (NRLS  "17). 

Same  as  NRLS 
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Symbols 

Definition 

Units 

NSIT 

Number  of  energy  levels  in  the  source 
spectrum* 

NSLS 

Number  of  energy  levels  in  the  source 
spectrum  of  the  transport  table*  (NSLS 
-  13) 

NSQR 

Number  of  source  altitudes* 

NSSI 

Same  as  NSLS 

NSST 

Option  to  specify  input  source  spectrum 
different  from  transport  table* 

NTRS 

Same  as  NR 

NUTS 

Checks  to  see  if  output  and  input  spectra 
have  coincident  energy  boundaries. 

ONNTR(l) 

The  neutron  flux  in  the  I  energy 
group  of  the  output  spectrum 

nt/fta 

OSRID(I) 

The  upper  boundary  of  the  Ith  output 
energy  group 

Mev 

PCTNT(I) 

fraction  of  neutrons  in  the  I*1*1  energy 
group  of  the  source  spectrum 

PHI 

Given  one-Msv  equivalent  flux  value ,  it 
is  desired  to  find  the  range  at  which  this 
flux  occurs. 

• 

PNTR(l) 

The  number  of  neutrons  in  the  Ith  energy 
group  of  the  output  spectrum* 

PSI(X,  I,  J) 

The  fraction  of  neutrons  scattered  frosLthe 
l'in  energy  group  at  the  source  to  the  «rn 
energy  group  at  the  receptor,  which  are 
separated  by  a  distance  RGD(K).  (This  infor¬ 
mation  is  Stored  in  th9  transport  tables). 

RADS(l) 

Response  factors,  nt/cm*  to  Rads  for  the 
energy  group 

RADS/(nt/cma) 

RDEN 

Density  of  air  at  the  receptor 

lbs/ft3 

RQD(I) 

iL 

The  I  distance  in  the  transport  table 

ft 

SDEN 

Density  of  air  at  the  source 

lbs/ft3 
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aynbols 

Definition 

Units 

SDNR 

(*  5  <jr  from  source  to  receptor 

lbs/ft1 

SLora(i) 

Used  to  calculate  direct  flux*  It  is 
the  of  lo8  ^fl'ux^  rang®  for 
the  I  energjr  group. 

1/ft 

SNKT 

Nuafcer  of  kilotons  being  investigated 

kT 

SNNTR(I) 

th 

Neutron  flux  at  recejfcoar  in  the  I 
energy  group  of  the  transport  table 
spectrun* 

nt/ft* 

SRID(I) 

Ith  value  of  transport  table  recejtttr 
energy  groups 

Kerr 

SSID(I) 

jth  Qf  transport  table  source 

energy  groups 

Kerr 

sspctr(i) 

Ith  value  of  input  source  energy 
groups 

Kerr 

tanntr(m,n) 

Intermediate  flux  calculation. 

TNNTS 

The  total  number  of  source  neutrons* 

26 


: Print  Air  Trans¬ 
port  Table? 

!YES _ NO 


FLOW  CHART 


Compute  by  log-linear  inter¬ 
polation  flux/nt  for  equr 


Print  j 
Table  j 


Read  Input 
Parameters 


Is  Source  Spectrum 
different  than  Air 
Transport  Spectrum  ? 
YES  NO 


Convert  source 
spectrum  to  air 
transport  spec¬ 
trum 


Multiply  by- 
total  spectrum 


Option  to  print 
transport  spectrum 
for  equr  ? 

YES  NO 


Print  | 


Is  output  spec¬ 
trum  different 
than  air  trans¬ 
port  spectrum  ? 
YES  NO 


Convert  transport 
spectrum  to  output 
spectrum 


Select  source 
altitude 


Option  ? 
;yes  NO 


